Phase Diagrams of Spinor Bose Gases 
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Using effective field theories dictated by the symmetry of the system, as well as microscopic 
considerations, we map out the magnetic coupling-temperature phase diagrams of spin-1 Bose gases 
in both two- and three-dimensions. We also determine the nature of all phase boundaries, and 
critical properties in the case of 2nd order phase transitions at both zero and finite temperatures. 



Spinor Bose condensates have been of strong interest 
for over a Due to the presence of hy- 

perfine degrees of freedom, such condensates have very 
rich internal structures and support a variety of novel 
topological defect s 0, IE El- In particular, depending on 
the values of two-particle s-wave scattering lengths in the 
singlet (S = 0) and quintuplet (S = 2) channels, ao and 
ci2, the system can form either nematic (polar) or ferro- 
magnetic spinor condensates^. In the meantime effects 
of thermal fluctuations are also much stronger in these 
systems, and it was shown recently[7] that the nematic 
(or polar) spinor condensate of spin-1 atoms is unstable 
at any finite temperature in two-dimension (2D), yield- 
ing to a pair condensate with power-law superfluid order 
until a Kosterlitz-Thouless (KT) transition temperature 
is reached. "Pairing" here is due to the thermal disor- 
dering of the nematic order parameter, instead of forma- 
tion of bound states between atoms, and the resultant 
state is found to have topological spin order Q. Above 
the KT temperature the system loses superfluidity alto- 
gether. Related workQ has focused on the case with 
easy-plane anisotropy, also in 2D, and found a multitude 
of KT transitions. Comparatively speaking, there has 
been much less discussion on the fate of ferromagnetic 
condensates, as well as what happens in three-dimension 
(3D). 



In the present work we clarify the phase diagrams of 
spin-1 Bose gases with spin-rotation symmetry, in both 
2D and 3D, and treat the nematic and ferromagnetic con- 
densates on equal footing. Our results are summarized in 
two phase diagrams, Fig. 1 and Fig. 2, for the 2D and 3D 
cases respectively. Of crucial importance in our consider- 
ation are the existence of a high-symmetry boundary that 
separates the nematic and ferromagnetic phases, and a 
quantum critical point separating the nematic conden- 
sate state from a pair condensate state. The rest of the 
paper will be devoted to elaborating on the physics asso- 
ciated with the various phases and phase transitions in 
these two phase diagrams. We will also briefly comment 
on what happens in one-dimension (ID) toward the end. 



The microscopic Hamiltonian describing spin-1 atoms 



with short-range interaction takes the formQ: 

/- ^2 



(1) 



where ip is the boson field, M is atom mass, a, a', 6, b' 
are spin indices running from —1 to +1, the vector F has 
three components which are the 3 spin-1 matrices, and 
repeated indices are summed over, cq and C2 are two- 
atom s-wave interaction constants; specifically, cq — 2c2 
and Co + C2 correspond to the interaction strength in the 
total spin and 2 channels respectively Stability of 
the system requires Co > 0, and cq > |c2 1 when C2 < 0, 
although the actual value of cq is not very important in 
determining which phase the system is in; we thus sim- 
ply assume it takes some sufficiently large value such that 
the two-body interaction is always repulsive in the quin- 
tuplet channel (thus 02 > 0) and do not discuss it further. 
We have also neglected the trap potential for simplicity, 
and assume instead the system is translationally invari- 
ant and in the thermodynamic limit. The Hamiltonian 
(HJ) possesses an U(l)xSU(2) symmetry, where U(l) cor- 
responds to charge conservation, and SU(2) corresponds 
to spin-rotation symmetry. 

The case with C2 = is very special. In this case the 
interaction, and therefore the entire Hamiltonian ([T]) is 
spin-independent. As a result the Hamiltonian JT]) has a 
much higher U(l)xSU(3) symmetry (llj. which one can 
understand by treating the 3 internal spin states as the 
fundamental representation of the SU(3) group. This 
case is of special importance in our discussion below. Ex- 
perimentally, this corresponds to the case that the scat- 
tering lengthes in the singlet and quintuplet channels arc 
the same: ao = 02 > 0. 

The spin structure of the spinor condensate is deter- 
mined by C2 [3j . For C2 < it is energetically favorable to 
have all atom spins aligned ferromagnetically; the ground 
state is a fully magnetized ferromagnet. For C2 > 0, spin 
polarization is suppressed and nematic spin states be- 
come energetically favorable. In particular, as we show 
below, when C2 becomes big enough, singlet bound states 
or molecules may form, and spin order disappears. 

Phases and Quantum Phase Transitions at Zero Tem- 
perature - We start by considering the phases and phase 
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FIG. 1: (Color online) Phase diagram of a spin-1 Bose gas in 
2D. At T = there are three phases: Ferromagnetic Conden- 
sate, Nematic Condensate and Pair Condensate, separated by 
quantum phase transitions at C2 = and C2 = C2 respectively. 
At finite T only the Pair Condensate phase survives, up to a 
Kosterlitz-Thouless transition, above which the system is in 
a normal phase with no broken symmetry. 



transitions at T — 0, where our analysis applies to both 
3D and 2D. We expect 3 different phases. 

(i) Ferromagnetic condensate (FMC) phase with C2 < 
0, where < a 2 < arj. In this phase atoms condense into 
a spin-polarized state: (fmc = e l6 (l, 0, 0) T , up to an 
arbitrary global rotation and gauge transformation. The 
order parameter is parameterized by (n, e l6 ), where 9 is 
the superfluid phase angle and n is a unit vector indi- 
cating the direction of spin polarization; thus the FMC 
order parameter lives in manifold 



M FMC = U(l) x S 2 



(2) 



Here S 2 (or 2-sphere) is the space of n, U(l) (also referred 
to as S 1 in the literature) is the space of 9. The dimen- 
sionality of this order parameter space is thus 2+1 = 3. 
The ground state spin quantum number is Stot = N, 
where N is the atom number. So the ground state is 
fully magnetized. 

(ii) Nematic condensate (NMC) phase with < C2 < 
cJj, where 0,2 > &o > — 00. In this phase atoms con- 
dense into a nematic state: (nmc = e (0, 1,0) T , up to 
an arbitrary global rotation and gauge transformation. 
The order parameter can be parameterized as a product 
ne %e , where 9 is the superfluid phase angle and n is a 
unit vector indicating the direction of spin nematic. It 
is extremely important to identify (9, n) and (9 + ir, — n) 
as the same point in order parameter space [5|, which is 
implied in our parameterization, ne ie . As a result the 
nematic spinor condensate order parameter lives in the 
coset manifold space 0,13] 



where the Z 2 reduction corresponds to the identification 
of (9, n) with (9+ir, — n). The dimensionality of the order 
parameter is still 3. 

It is clear that the boundary separating the FMC and 
NMC phases is at C2 = 0, where ao = a 2 . This is a 1st or- 
der phase boundary across which magnetization changes 
by N. 

(iii) As C2 increases, the two-particle interaction turns 
attractive in the singlet channel (and ao turns negative) , 
and beyond certain threshold bound states or singlet 
molecules can form (where ao turns back to be posi- 
tive). In this case the ground state is a Bose conden- 
sate of singlet molecules, which has no spin structure [lij]. 
This molecular, or pair condensate (PC) state breaks the 
U(l) symmetry associate with charge conservation, but 
respects the SU(2) spin rotation symmetry; its order pa- 
rameter is e 10 , where 9' is the superfluid phase angle of 
the molecules. Thus the PC order parameter space is 
simply 



M, 



PC 



mi) 



(4) 



which is 1-dimensional. As a result there must be a quan- 
tum phase transition between the NMC state and the PC 
state at T = 0, driven by increasing quantum fluctuation 
parameterized by C2. We identify the QPT point as c|, 
which turns out to be a quantum critical point (QCP), 
as we show below. In the dilute limit, c 2 corresponds 
to the Feshbach resonance in the singlet channel where 
\a \ -> 00. 

To study the critical properties of this QCP, we write 
down the quantum (imaginary time) effective action that 
describes the NMC and PC phases on equal footing: 
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(5) 



M NMC = [U(l) x S 2 ]/Z 2l 



(3) 



where the first two terms is the action for the NMC[f|, 
the 3rd term is the action for the PC, and last term 
describes the conversion between two atoms and a sin- 
glet pair. Such "two-channel" models have been used 
to study pairing induced QPT in spinless[l3| and spin- 
1/2[1J| bosons. The form of this effective action is essen- 
tially determined by the symmetry properties of these 
two phases, and their order parameter dynamics or col- 
lective mode spectra. We have neglected terms consistent 
with symmetry but involve more gradients, which are ir- 
relevant in the long-distance and low- frequency limit. 

It should be noted that the PC order parameter 9' is 
ordered in both the nematic and pair condensates states; 



3 



it thus has a non-zero expectation value at the tran- 
sition, with fluctuations that can be safely integrated 
out 3, 14 1 . Thus for the description of this transition, 
we can choose (without loss of generality) 9' — 0, and the 
last term of Eq. ([5]) reduces to J' J drdr cos(20). This 
reduces 9 to an Ising or Z 2 -like variable, as the expression 
above takes minima at 9 — and 9 = 7r[l^|. However be- 
cause (9, n) and (9 + ir, — n) correspond to the same point 
in the NMC order parameter space, we have the freedom 
of making the "gauge-fixing" choice of 9 = 0; this in turn 
removes the redundancy in the space U(l) x S 2 , thus n 
now lives on the full S 2 space after this "gauge-fixing" 
procedure. 

After integrating out massive fluctuations of 9, we ar- 
rive at the appropriate field theory for the QPT between 
NMC and PC phases: 



o _ Pn 
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(6) 



where n is a unit vector of 3D spin space that lives in the 
full S 2 space. This is nothing but the 0(3) non- linear a 
model (NLcrM). We thus conclude the universality class 
of the QPT at c 2 is that of d + 1-dimensional 0(3) or 
Heisenberg transition where d = 2 or 3 is the spatial 
dimension of the system. 

We note that this result is quite unusual, and the pres- 
ence of superfluidity on both sides of the transition is 
of crucial importance. Normally one expects a 1st or- 
der transition between nematic and isotropic phases [la] , 
when no superfluidity is involved. In insulating quantum 
magnets a continuous transition is possible, but it is of a 
more exotic type [17 1. 

In the theory ([5]) we have neglected the atom-molecule 
and molecule-molecule interactions. In the case of spin- 



less bosonic atoms [13j such interactions are generally at- 
tractive, which render the PC phase unstable in that case. 
In the following we show this is not the case here. This is 
most easily demonstrated by going to the strong pairing 
limit c 2 — > oo, and use a lattice regularization in which 
the integrals in Eq. (|TJ) are turned into summations over 
lattice sites. In this case the interaction terms in (JTJ) take 
the form 

V = [(co/2)n i (n i - 1) + (c 2 /2)(s 4 2 - 2n*)], (7) 

i 

where hi and Si are the total atom number and spin 
operators of site i respectively. In the strong pairing 
limit the molecule is simply a singlet pair of two atoms 
on a given site, whose energy is E pa i r = cq — 2c2, while 
a single atom has energy E at0 m — 0. Now consider the 
interaction energy of a molecule and an atom occupying 
the same site (or equivalently, 3 atoms on one site with 



total spin s = 1): E 3 = 3c — 2c 2 > E, 



pair 



E a t 



Similarly, consider two molecules on the same site (or 4 
atoms forming a total singlet): £4 = 6cq — 4c 2 > 2E pa i r . 
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FIG. 2: (Color online) Phase diagram of a spin-1 Bose gas 
in 3D. Thick solid lines are 1st order phase boundaries, while 
dashed lines are 2nd order phase boundaries, with universality 
classes indicated with the same color. 



We thus conclude that the atom-molecule and molecule- 
molecule interactions are both repulsive in the present 
case, and thus the PC phase is stable. 

Finite Temperature Phases and Phase Transitions in 
2D - It is well-known that thermal fluctuations sup- 
presses ferromagnetic or nematic spin order in 2D[16|. 
As a result among the three phases at T = 0, only the 
PC phase survives at finite T (with power-law long range 
order), which eventually yields to the normal phase via 
a KT transition. All these are reflected in the phase di- 
agram, Fig. 1. It is interesting to note that within the 
PC phase and at very low T, as one goes from small to 
large c 2 , there is a crossover similar to the BCS-BEC 
crossover in spin- 1/2 fermionic superfluids, in the follow- 
ing sense. In this phase one can identify the pair size £ 
as correlation length of the spin nematic order parame- 
ter. At small c 2 we have £ 3> I, where I is the average 
interatomic distance, very much like what happens in a 
weak-coupling BCS fermionic superfluid. For large c 2 
atoms form closely bound singlet molecules, thus £ is the 
size of the molecule, and we have £ <C I, similar to the 
BEC regime of the fermionic superfluid. The crossover 
happens near c 2 . 

Finite Temperature Phases and Phase Transitions in 
3D - The 3D phase diagram is considerably richer, as all 
three phases at T — are stable at low-T. Similar to 
the 2D case, the PC phase yields to the normal phase 
via a 3D XY (or 0(2)) transition as T increases. For 



c 2 



< 



(in the NMC phase but near the QCP), we expect 
the nematic order parameter to be suppressed and the 
NMC phase first yields to the PC and then to normal 
phase as T increases. Using analysis similar to the QPT 
between NMC and PC phases, we conclude the classical 
NMC to PC transition is in the universality of 3D 0(3) 
transition. The transition temperature T c near the QCP 
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is determined by the critical properties of the QCP[l8j: 

T D <x {&-(*)'", (8) 

where in the present case the dynamical exponent z = 1 , 
and because d+ 1 = 4 is the upper critical dimension, the 
QCP correlation length exponent v — 1/2 but there may 
be logarithmic correction to ((8|). For reasons to be elab- 
orated below, we expect the 0(3) and XY phase bound- 
aries to cross at c' 2 with < c' 2 < c|, so for < c 2 < c' 2 , 
we expect the NMC phases to first yield to a normal 
phase with nematic spin order via an XY transition, 
and then to the isotropic normal phase. The nematic 
to isotropic transition is expected to be 1st order 

On the ferromagnetic side (c 2 < 0), condensates (in- 
cluding multiparticle condensates if present) must also 
have ferromagnetic spin order. Thus the order of transi- 
tions must be that the FMC phase first yields to a normal 
ferromagnetic phase via an XY transition, and then to 
the isotropic normal phase via an 0(3) transition. This 
persists all the way to the SU(3) line of c 2 = 0. The 
separation of magnetic and superfluid transition temper- 
atures can be easily established by a perturbative treat- 
ment of Co and c 2 . 

Due to the enhanced symmetry, along the c 2 = line 
the order parameter lives in a much larger space, which 
is a 5-dimensional manifold: 

M C2=Q = U(l) x {SU(3)/[U(1) x SU(2)}}, (9) 

where the first factor U(l) is the same as in previous 
cases, while the 2nd factor represents the 3-component 
spinor order parameter <j> with the constraint cfvcj) = 1 and 
excluding its overall phase; the SU(2) reduction there is 
the subgroup of SU(3) that leaves <f> invariant. 

For the same reason as in the c 2 < case, as T in- 
creases, we expect the U(l) symmetry to be restored first, 
and then the SU(3) symmetry. The effective action for 
the SU(3) critical point is the CP2 model: 

S[4>] = ^J dr[(V + iA)flt • [(V + iA)fl. (10) 

The high symmetry also guarantees the c 2 = line is 
a 1st order boundary separating ferromagnetic and ne- 
matic phases, with or without condensates. Furthermore, 
the XY transition line and isotropic to ferromagnetc or 
nematic transition line must be continuous at c 2 = 0. As 
a result for small positive c 2 there must be a nematic 
normal phase separating the NMC and isotropic normal 
phases. All these are summarized in the 3D phase dia- 
gram, Fig. 2. 

Phases in ID - Due to strong thermal fluctuations, 
the only stable phase at finite T is the isotropic normal 
phase. At T = the NMC phase is destroyed by quantum 
fluctuations, thus only the FMC and PC phases remain, 
separated by the 1st order QPT at c 2 = 0. The ferro- 
magnetic order parameter is a conserved quantity and 



thus not affected by the enhanced quantum fluctation, 
but the superfluid order parameters only have power-law 
long-range order in both phases. 
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